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ABSTRACT 

The active control of acoustic pressure in a 2-D cavity with a flexible boundary (a beam) 
is considered. Specifically, this control is implemented via piezoceramic patches on the beam 
which produce pure bending moments. The incorporation of the feedback control in this 
manner leads to a system with an unbounded input term. Approximation methods in the 
context of an LQR state space formulation are discussed and numerical results demonstrating 
the effectiveness of this approach in computing feedback controls for noise reduction are 
presented. 
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1 Introduction 


In recent years, the development of new fuel efficient turboprop engines has motivated the 
development of a comprehensive active control methodology for interior pressure field cham- 
bers. The active control of noise in this setting has been studied both in a frequency domain 
setting [14, 19] and from an infinite dimensional state space time domain approach (PDE 
approach) [2, 6, 7, 12] with techniques often centering around the generation of an appropri- 
ate secondary pressure wave which optimally interferes with the offending primary pressure 
wave. Here however, we consider a time domain state space formulation in which the active 
control is implemented via piezoceramic patches which are imbedded in the boundary of the 
acoustic cavity. 

The example we consider consists of an exterior noise source which is separated from 
an interior chamber by an active wall or plate. This plate transmits noise or vibrations 
from the exterior field to the interior cavity via fluid/structure interactions thus leading to 
the formulation of a system of partial differential equations consisting of an acoustic wave 
equation coupled with elasticity equations for the plate. The control is implemented in the 
example via piezoceramic patches on the plate which are excited in a manner so as to produce 
pure bending moments. It should be noted that the incorporation of the feedback control 
in this manner leads to a system with an unbounded input term. Experiments are being 
designed and carried out at NASA Langley Research Center in which the interior cavity is 
taken to be cylindrical with a circular active plate and sectorial patches. 

As a first step toward developing an effective linear quadratic regulator (LQR) state space 
control methodology for near field acoustic problems of this type, it is useful to consider a 
simplified but typical model consisting of a 2-D interior cavity with an active beam at one 
end (see Figure 1). Here T represents a perturbing force on the beam due to an exterior 
noise source. This in turn causes fluctuations in the interior acoustic pressure field and hence 
unwanted noise. The goal in the control problem is to optimally reduce the interior pressure 
deviations by effecting a force distribution on the beam that decouples the cavity acoustic 
response. 

In Section 2, a model set of differential equations for the problem is given and the math- 
ematical framework needed to pose the control system in an abstract Cauchy formulation is 
presented. Section 3 contains a brief discussion of the theory of finite and infinite dimensional 
periodic optimal control problems while Section 4 is devoted to the general finite dimensional 
approximation of the control problem. Specific approximation schemes are discussed in the 
fifth section and examples demonstrating the viability of the method are presented. 
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Figure 1 . Acoustic chamber with piezoceramic patches. 


2 Mathematical Model 

When describing acoustic wave motion in a fluid, it is useful to introduce a velocity potential 
<f> which is a complex- valued function satisfying v(t,x,y ) = —V<f>(t,x,y) where v denotes the 
fluid’s velocity [ 15 , 16 ]. If the equilibrium density of the fluid is given by />/, the acoustic 
pressure p (the deviation from the mean pressure at equilibrium) is related to this velocity 
potential by p(t,x,y) = pf<j) t (t,x,y). For acoustic waves with small amplitude, both the 
potential and the pressure satisfy the undamped first order wave equation with uniform 
speed of sound c in the fluid; hence 

<f>u = c 2 A<f> (x,y)€fi(f),*>0. 

The boundaries on three sides of the variable cavity tt(t) are taken to be “hard” walls thus 
leading to the zero normal velocity boundary conditions 

V<£-n = 0 (x,y)er,f>0 
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where n is the outer normal. It is assumed that the perturbable boundary consists of an 
impenetrable fixed-end Euler- Bernoulli beam with Kelvin- Voigt damping. If w(t,x) is used 
to denote the transverse displacement of the beam with linear mass density pb , the equations 
of motion are 


Q 2 0 < x < a , 

pbWtt + —M{t,x) = -pj<t>t(^X,w(t,x)) + f{t,x) t> Q , 

w(i,0) = = w(t,a) = = 0 t> 0 , 


( 2 . 1 ) 


where M(t,x ) is the internal moment and / is the external applied force due to pressure from 
the exterior noise field. For an uncontrolled beam with Kelvin- Voigt damping, the moment 
contains both strain and strain rate components and is given by 


M(t,x) = El 


d 2 w 

Je 2 


+ CpI 


d 3 w 

dx 2 dt 


The final coupling equation is the continuity of velocity condition 

w t {t,x) = V<f>(t,x,w(t,x))-h , 0<x<a,t>0 (2.2) 

which results from the assumption that the beam is impenetrable to fluid. Under an as- 
sumption of small displacements {w(t,x) = w{t,x) + 8 where w = 0) which is inherent in 
the Euler-Bernoulli formulation, the beam equation in (2.1) can be approximated by 


r\2 

PbWtt + ■— M{t,x ) = -p/fot(i,z, 0) + cf>ty{t,x,Q)w] + f(t,x) 

while (2.2) can be approximated by 

Wt(t , x) = V^>(U x, 0) • n + (V <j> y (t , x, 0)u>) • h . 

To first order, these last two equations can be approximated by dropping the higher order 
terms -/>/&»(*, x, 0)u> and (V^(t,*,0)u>) • n. Then upon approximating the domain Sl(t) 
by the fixed domain Q = [0, a] x [0,^], we obtain the approximate uncontrolled model 


(j) tt — c 2 A^> (x, y) € 0 ,< > 0 , 

V^> • n = 0 (x,y) € T ,t > 0 , 

^(t,x,0) = -wt{t,x) 0<x<a,t>0, 

oy 

pbwu + (■ ei I ? + *>*£&) = + f{t,x) 

w (t, 0) = ^(t, 0) = w(t, a) = ^(t, a) = 0 t > 0 , 

(j>{ 0, x, y ) = <t>o{x, y) , w(0, * ) = w o ( x ) 

<f> t (0,x,y) = <f>i{x,y) , «h(0,x) = u>i(x) . 


0 < x < a , 
t > 0 , 


(2.3) 
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For control of structural vibrations and the acoustic pressure field in this model, s piezo- 
ceramic patches are attached to the beam as shown in Figure 1. These patches are excited 
in a manner so as to produce pure bending moments ([8, 9, 11]) (see Figure 2). If H is used 
to denote the Heaviside function, the model for the controlled beam can be written as 

, d 2 ( d 2 w T d 3 w \ 

PbWtt + w[ E 1 W +CDl ^t) +P ^ X ^ 

d 2 / K B k s \ ( 2 ‘ 4 ) 

_ ~dx 2 T u «(0 \H{, X ~ a »i) ~ H(x — a, -a )] ) + . 

Here n,(t) is the voltage applied to the i tfi patch, K B is a parameter which depends on the 
geometry and piezoceramic material properties, T is the patch thickness and A; is a material 
constant (see [8, 9]). It should be noted that the incorporation of (2.4) into (2.3) leads to a 
system with an unbounded input term since it involves the second derivative of the Heaviside 
function. 


+ 



Figure 2. Piezoceramic patch excitation. 


To formulate this problem in the context of existing infinite dimensional control theoretic 
results, it is advantageous to pose the control system in an abstract Cauchy formulation. To 
accomplish this, the state is taken to be 2 = ( <f>,w ) in the Hilbert space H = Z 2 (Q) x T 2 (r 0 ) 
with the energy inner product 



Here L 2 (Q) is the quotient space of L 2 over the constant functions. We also define the 
Hilbert space V = H 1 (Q) x //^(To) where /A 1 (H) is the quotient space of II 1 over the 
constant functions and // 0 2 (r 0 ) = {V> € H 2 (T 0 ) : ^>(x) = = 0 at ar = 0,a}. The V inner 

product is taken as (here and below we use the notation D — ^-1 

OX ' 


4 



Following the ideas used in the theoretical results in [3, 4], we consider the Gelfand triple 
V <_► H V* with pivot space H and define sesquilinear forms <r, : V x V -> <C , i = 1,2 

by . 

<7^$, ^)= / p } V<j>-V(duj + / EID 2 wD 2 r)d r y , 

Jo Jr 0 

cr 2 ($,'I') = [ {c d ID 2 wD 2 t] + p f {fa-wO}d 7 

where $ = and 'f r = (£,* 7 ) are in V. It can be easily argued that the sesquilinear 

forms satisfy the continuity and coercivity conditions 

Re > Ci|$|v , 

<C2|*lv|*k » 

Re <r 2 ($, $) > c 3 ( D 2 w , D 2w ) L 2 (ro) = c 3l w l^(r 0 ) > 

<C4|*M*|V 

(for detailed arguments in a similar setting, see [1]). The control operator B € C(U,V ) is 

defined by g 

(B«, «> v . v = [ El ^ u i( H « - 

•To 1 i = 1 

for where fl-y(a) = i/(x - oq),* = 1,2, •••,«, j = 1,2 and (vJv.v 1S the usual 

duality pairing. , , ,, 

Finally, for F = (0, f/pb) we can write the control system in weak or variational torm 

( z tt (t ), $) v .,v + + cr i( z (^)> 1 ®) = ( Bu(t ) + F, y) v *,v ^-5) 

for $ in V. The state is given by z{t) = W, •, •), w(t, •)) in V ^ //. Since cr x and <r 2 are 
bounded, we can define operators A X ,A 2 G C(V : V*) by 

(A^^) v . y = <r i (^,^) 

for i = 1,2. This then yields the system 

z u {t ) + A 2 z t {t) + A x z(t) = Bu{t) + F 


i n y*. ....... 

Continuing with our abstract formulation, we next write the system in first order torm. 
To accomplish this, define the product spaces V = V xV and H = V x H with the norms 


= |4>| 2 V + |*& 


and 


|(*,tf)lv = l*lv + l*lv- 
For x — ($, and © = (T, A), the sesquilinear form cr : V x V — > © is then defined by 

cr((T, A), ($, tf)) = - (A, $) v + cr x (T, $) + cr 2 ( A, '5) . (2.6) 
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Since the duality product {•, -) v .y is the unique extension by continuity of the scalar product 
{•> ' )h f rom H x V to F* x V, it follows that for appropriate restrictions on 0 we can write 

a(0,x) = <r((T, A), ($,$)) = -(A,*) v + (^ 1 T,$) K . tV + (A 2 A,«) v . >v , 

= -(A.^y + ^T + ^A,^ 

= ((-A 1 A l T + A 2 A), (*,*))„ 

= (-^0,x) w - 

The operator A : 'H — » T~i is given by 



where dom A = {0 = (T, A) € H : A € V t Ai T + A 2 A € i/}, Ai and A 2 are the operators 
defined by a i and <r 2 , respectively, and the above calculations hold for 0 £ dom A (see [1] 
for further examples concerning the definitions of operators and domains in this manner). 

To write the first order system in weak or variational form, let Z(t) — (z{t), z t (t)), 
^"(0 = ((W)), an( ^ Bu(t) = (0, Bu(t)). The weak form of the system is then 

{■£<(*)> x)v,v + X) = (Bu(t) + T(t), x)v,v (2-8) 

for x € V. Formally, this is equivalent to the system 

Zt(t) = AZ(t) + Bu(t) + F(t) ( 2 . 9 ) 

in 7 i where A is given in (2.7). 


3 Periodic Control Problems 

As noted in the introduction, our control problem is motivated by the desire to reduce 
cavity pressure fluctuations resulting from the perturbing noise J~ . In many applications, 
it is reasonable to assume that J~ is periodic with period r; hence an important problem 
of interest (e.g., see [6]) for the system (2.9) is an LQR problem for a periodic disturbing 
force T. This can be formulated as the problem of finding u € T 2 (0,t; U) which minimizes 
a quadratic cost functional of the form 

J(U) = llo Z ^ n + <*“(*). dt 

subject to (2.9) with 2^(0) = Z(t). Since Z = (<f>, w, (f> t , w t ) T , the operator Q can be chosen 
so as to emphasize the minimization of particular state variables as well as to create windows 
that can be used to decrease state variations of certain frequencies. The control space U is 
taken to be HT if s patches are used in the model, and it is assumed that the operator R is 
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an s x -s diagonal matrix where r,-,- > 0, i = 1, ■ ■ • , s is the weight on the controlling voltage 
into the i th patch. In the case that B is bounded on Ti, a complete feedback theory for 
this problem can be given as discussed in [10]. Under usual stabilizability and detectability 
assumptions on the system as well as standard assumptions on Q , the optimal control is 
given by 

u(t) = -R~ x B m [HZ{t) - r(t)} 

where II is the unique nonnegative self-adjoint solution of the algebraic Riccati equation 

A*n + nA-nBR- 1 B*n + Q = o . ( 3 . 1 ) 

Here r is the unique r-periodic solution of 

r(t) + ( A * - n BR~ l B*)r{t) - n T{t) = 0 (3.2) 

and the optimal trajectory Z is the solution of 

Z{t) = {A- BR~ 1 B*Tl)Z(t) + BR~ l B*r{t) + T{t) . 

These equations (in particular (3.1), (3.2)) are infinite dimensional (i.e., in Ti) and hence 
approximation techniques are required to obtain approximate feedback gains. Using a stan- 
dard Galerkin approach, one typically chooses a sequence of finite dimensional subspaces 
T~L N C 77 with projections V N : H — > H N and defines an approximating problem in 'H N of 
minimizing 

J N (u) = 1 £ {(Q N Z N {t),Z N (t)) n + ( Ru(t ), <1(()V} it 
subject to an approximating system 

Z N (t) = A N Z N (t) + B N u(t) + T N {t) 

Z N ( 0) = Z n (t) = V N Z( 0). 

The solutions are given by 

u N {t ) = -R-'B^^ZV^t) - r N (t )] 

Z N (t) = (A* - B N R~ l B N *n N )Z N (t) + B N R~ 1 B N *r N (t) + Z N {t) 
where is the unique nonnegative self-adjoint solution of 

. 4 ** 11 * + n n a n - u n b n r- 1 b n *u n + q n = o 


and r* is the unique r-periodic solution of 

r N (t ) + (^** - n N B N R~ 1 B N *)r N (t) - n N f N (t) = 0. 

In order to guarantee the convergence U N 'P N Z —* HZ for Z €. W, r N (t ) — > r(t ), and 
hence the convergence of u N (t) to u(t), it is sufficient to impose various conditions on the 
original and approximation systems. These hypotheses include convergence requirements for 
the uncontrolled problem as well as the requirement that the approximation systems preserve 
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stabilizability and detectability margins uniformly. A fully developed theory (see [ 4 ]) is 
available for the case that T = 0 (in this case the tracking variable r does not appear in the 
solution) even in the case that B is unbounded in the sense formulated in Section 1, The 
theory in [ 4 ] requires rather strong damping assumptions on the second order system ( 2 . 5 ) 
in order to be applicable. Under appropriate assumptions, the techniques and ideas of [ 3 ] 
and [4] can be used to treat the case for T ^ 0 in both identification and feedback control 
problems. 


4 Finite Dimensional Approximation 

An advantageous feature of the state space approach for feedback control is that the optimal 
control can be implemented using various approximation techniques. To illustrate the ideas 
involved, let {B"}"_r i 1 denote the 1-D basis functions which are used to discretize the beam 
and let {B™}™^ m = (m x + 1) • (m y + 1) — 1, denote the 2 -D basis functions which are 
used in the cavity. The n — 1 and m dimensional approximating subspaces are then taken 
to be H£ = span and H™ = span {B™}™ x , respectively. Defining N = m + n — 1, 

the approximating state space is H N = H™ x and the product space for the first order 
system is 7 i N = H N x H N . The finite-dimensional approximation is then determined by 
restricting a to x where a is given in (2.6). This yields the operator A N : H N — > 'H N 
where 

I A n = 

and and are obtained by restricting o\ and cr 2 to II N x H N . We observe that the 
restriction of the infinite dimensional system ( 2 . 5 ) to the space H N x H N yields for $ = (£, 77) 

; + ».(*"(<). #) 

, = / / fld -1 . 

-To * *To 

When ^ is chosen in H N and the approximate beam and cavity solutions are taken to be 

f «,"((,*) = ixwBh*) 

: 1=1 

and 

m 

»=1 
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respectively, this yields the system 

M N y N (t) = A N y N (t) + B N u(t) + F N (t ) 
M N y N ( 0) = 

where 


V ->"(<) ) 


Here ^(t) = (<j>? (t), w? (t), ■ ■ ■ , w^_ l (t)) T denotes the N x 1 — 

(m + n — 1) x 1 approximate state vector coefficients while u(<) = (ui(t), • • • i/ s (f)) contains 
the 5 control variables. The full system has the form 


" M? 

0 

' 0*(O ' 

0 

M? _ 

. ^(0 . 

' 

0 

' 0 ^( 0 ) ' 

0 

m" _ 

_ 0*(O) _ 


r o m* l 

' 1 9 N (t) ' 


0 


0 

i 

i 

-‘5- 

1 

1 


+ 

b n 

) -+■ 

F N (t) _ 


with 


9 ? 

9 ? 


= diag[M^, Mx 2 \ , 
M 2 n = diag[M£, M%) , 
Ai = diag^,^] , 

0 


and 




B n = 




6 jV 

^21 




0 


F*(t) = 


0 

L Wo J 
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The component matrices are given by 

K] ( t = j a VB? ■ VB?dw , [M«] p . = £ D 2 B” D 2 B”d~y , 

KL, = i ^BTBTdu, , [Mg] p . = / pt pbB™B*d.~i , 

K) , t = l PI^B? ■ VB,”*., , [AS] R . = EID 2 BlD 2 B;d-t , 
K,L = - f PjB?BTd-r , [A«] = / p,B?B;i-, , 

‘'lo P } 1 o 

[< 2 ] pi = ^ c d ID 2 B? D 2 Bp d'y , 

[4"] . = P 2 EI^r~D 2 Bpd'y , 

1 j pj y ttil t y 

[^w] p = 

Moreover, the vectors <7^ = [<7n, <7(2] and = [</£!, <7^] have elements 


[jfi] , = J a v<h ■ VB?dw , [<&] p = / ro D 2 w 0 D 2 B;di , 
[?£], = jf , [s!*] p = l' w,b;^ . 


In all cases, the index ranges are k,t — 1, • ■ • , m and z,p = 1, • • • , n — L The patch index 
j ranges from 1 to s. It should be noted that the matrices 4^ and M N are symmetric and 
positive definite by construction. The matrix has a symmetric block and a skewsymmetric 
block and the eigenvalues of 4^ are real and nonnegative. 

With the bases {2?J l }" = r i 1 and chosen, the finite dimensional theory outlined in 

the last section holds with the various finite dimensional operators replaced by appropriate 
matrices. Specifically, the finite dimensional control problem is then to find u £ L 2 (0,r) 
which minimizes 

jN ( u ) = \j l . iV = m + n-l 

where Q N is nonnegative definite and y' solves 

»"(<) = dV(l) + B N u(t) + F N (l) 

N (4- 1 ) 

y { o) = Vo ■ 

Here A N = (M n )~ 1 A N , B N = [M N )~ l B N , F N (t) = (aP)' 1 F N (t) with the initial 
condition i)q — tjq . The optimal control is 

u N (t) = R- 1 (B n ) t [r w (0 - n N y N (t)] (4.2) 
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where n w is the solution to the algebraic Riccati equation 

( A*) T n" + n n a n - u n b n r-' ( b n ) t u n + q n = o . (4.3) 

Since Q N denotes the matrix representation for the operator Q N , a suitable choice for Q N is 

M? 0 


q n = v 


Mf 


where the diagonal matrix T> is given by 

V = diag [d 1 / m ,d 2 / n - 1 ,4/ m ,d 4 /"' 1 ] • 

Here I k , k = m, n — 1 , denotes a k x k identity and the parameters d, are chosen to enhance 
stability and performance of the feedback. The sxs diagonal matrix R contains the positive 
control weights and has entries ra,i = l,-”, s. For the regulator problem with periodic 
forcing function F^ (t), r N (t) solves the linear differential equation 


r N {t) = - [A n - B n R- 1 (B n ) t U n ] T r N (t) + U N F N (t) 
r N ( 0) = r N (r) 

while the optimal trajectory is the solution to the linear differential equation 

y N (t) = [A* - B N R-'(B N ) T n N ] y N (t) + B"R-'(B N ) T r N (t) + F N (t) 

y N ( 0 ) = y N {r) ■ 


(4.4) 


(4.5) 


5 Specific Approximations and Numerical Results 

We next turn to a discussion of specific choices of basis functions in the general formulation of 
the approximation schemes in the last section. We shall also present numerical results from 
related computations. When choosing bases for the finite dimensional subspaces H£ and 
H™ in a control setting, one must weigh criteria such as smoothness requirements, uniform 
preservation of exponential stability of approximating systems (see [5]), accuracy, sparsity 
of system matrices and ease of implementation. 

From energy considerations, it follows that the system (2.3) is dissipative; hence all the 
eigenvalues lie in the left half plane. The model of Section 1 includes no medium damping 
however, and hence the energy dissipation in the cavity results exclusively from the boundary 
(Kelvin- Voigt damping in the beam) thus making the system (2.3) only weakly damped. 
In spite of the lack of strong damping, numerical tests have indicated that when physically 
relevant parameters are used in the model, the system (2.3) is exponentially stable (the fixed- 
end boundary conditions on the beam make difficult a thorough analytical analysis of the 
eigenstructure). When considering various methods of discretizing the problem, one would 
like to choose schemes which uniformly preserve the exponential decay rate as the dimension 
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of the approximate system (4.1) increases. This can be easily checked by determining whether 
or not there exists a uniform margin for increasing N between the open loop eigenvalues of 
the system matrix A N in (4.1) and the imaginary axis. 

When considering the control problem, one is also concerned with the preservation of 
uniform stabilizability and detectability margins for the closed loop approximation systems. 
Hence care must also be taken so that approximation schemes are chosen so as to preserve 
a uniform margin between the closed loop eigenvalues of A N - B N R _1 (B N ) T U N and the 
imaginary axis. Numerical schemes which satisfy these various criteria will now be discussed. 

Cubic splines were used as a basis for H% since they satisfy the smoothness requirement 
as well as being easily implemented when adapting to the fixed-end boundary conditions 
and patch discretizations. For a given positive integer n, a uniform partition was taken with 

the gridpoints a:" = £a, i = 0, 1, ■ ■ • ,n. If is used to denote the standard cubic 

spline basis corresponding to this partition (see [18], page 79), then the basis functions for 
the beam discretization were taken to be 

£?" = Bq — 2B™ - 2B n _ x 

B™ — B* ; i = 2, 3, • • • , n - 2 

BU = K - 2 BU - 2 B: +1 

It is readily seen that these basis functions satisfy the essential boundary conditions; that is, 

£, n (0) = DB?(0 ) = H”(a) = DB?(a) = 0 

for i = 1, 2, ■ ■ • , n — 1. As mentioned previously, the corresponding n — 1 dimensional 
approximating subspace is then given by H% = span and the approximate beam 

solution is taken to be 

w N (t,x) = w?(t)B?( 2 r) . 

t=l 

With this choice of basis functions, the matrices and M ^ are easily constructed and 
are 7-banded. It should be noted that a Tau-Legendre discretization was also considered 
for the beam but had the disadvantage of the loss of four equations due to the constraints 
mandated by the fixed-end boundary conditions (see [13] for a discussion of Tau methods). 

The bases that were considered for the cavity discretization included tensored one- 
dimensional Legendre polynomials, tensored linear splines and finite elements. The methods 
of system formulation as well as the advantages and disadvantages of each can be summa- 
rized as follows. Consider first the Legendre basis. Let P?(x) and P-{y) denote the standard 
Legendre polynomials that have been scaled by transformation to the intervals [0, a] and 
[0,£], respectively. The basis functions for the cavity are then defined as 

Bij ( x ’>y) = P{ ( x )Pj{y ) f° r * = 0, 1 , • • • , m x , j — 0, 1, ■ ■ • , m y , i + j ^ 0, 

where m = ( m x + 1) • (m y + 1) — 1. The condition i -f j ^0 eliminates the constant function 
thus guaranteeing that the set of functions is suitable as a basis for the quotient space. For 
definiteness, the basis functions are ordered by assuming that i varies for each fixed j which 
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is analogous to a left to right, bottom to top ordering. Notice that because natural boundary 
conditions occur on all sides of the cavity, one does not have to employ a Tau method; that 
is, the method is simply a Galerkin scheme without modification of the basis elements to 
satisfy some essential boundary conditions. 

The component matrices Mg x and Mg x can then be succinctly described as follows. Let 
the fundamental (m r + 1) x ( m x +1) matrices M™ and K™ be defined as 

[AC] ,i= f WW* 

J o 

[/C] 0 -= r DP?(x)DP*(x)dx 

J 0 

with similar definitions for Jiff, I<g. Using the tensor properties of the 2-D basis, we can 
form the matrices Mg and Mg 1 defined by 

M([ = M? 0 I<T + K? ® 

Mg = %Mg® m™ . 

The ordering in the above definition depends on the ordering of the basis functions. The 
matrices Mg t and Mg are obtained by removing the first row and first column of Mg and 
Mg to reflect the deletion of the constant function from the basis set. Note that with 
this definition, both matrices are very easily constructed and that the mass matrix M 21 is 
diagonal; hence the inverse is trivial to calculate. Although the matrix Mg is not sparse, it 
has a well-defined structure due to its tensor product nature and the fact that Mg and Mg 
are diagonal. It too can be efficiently inverted when one takes advantage of this structure. In 
the case that pj is constant, the stiffness matrix Ag can be constructed in the same manner 
as Mg and the tensor product structure can be used advantageously both when solving the 
Riccati equation (4.3) and the ODE systems (4.4) and (4.5). 

In order to use a tensored linear spline or finite element basis in the cavity, some constraint 
must be applied in order to guarantee that Hg is a quotient space (one cannot simply drop 
the constant function as was done with the tensored Legendre polynomials). One such 
constraint which is commonly used is the requirement that 

/ 4> n ( t,x,y)du> = 0 . 

J n 

If is used to denote the standard tensor product linear spline basis (see page 129 

of [18]), then the integral constraint leads to the quotient basis {Bg }™ 2 where 

Bg{x,y) = B™(x,y) - 4a i R 1 m (x,j/) 

with aj = /j, 1 1 depending upon whether the function B™ is a corner basis function, a 

side basis function or an interior basis function, respectively. As a result of the modifica- 
tions needed to obtain a quotient space basis with the tensored linear splines, the matrices 
Mg, Mg and An are full and hence one loses the structural advantages obtained with the 
Legendre basis. 
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Similar modifications must be made when using a finite element basis in a quotient space 
with the result that the system matrices are also full in that case. Moreover, the lower 
order accuracy of the splines and finite elements necessitates the use of a larger number of 
basis functions and hence larger matrices in order to match the accuracy of the Legendre 
polynomials. The fact that the Legendre basis yields smaller, structured matrices than 
those obtained with the linear splines and finite elements is important but not crucial in the 
problem under consideration since the cavity is only two dimensional and hence matrix sizes 
are reasonably small. This issue will become much more critical when considering the 3-D 
problem of interest because of the large matrix sizes which will be encountered. 

As discussed earlier, a final item w r hich should be considered when choosing a means 
of discretizing the control problem is whether or not the approximation scheme effects a 
uniform preservation of exponential stability for the open and closed loop approximating 
systems. This issue is illustrated by the results in the Example 5.1. 

The problem under consideration in Examples 5.1 and 5.2 is 


(j>u = c 2 A4> (x, y) € D , t > 0 , 

Vcf>-h = 0 (x, y) £ T ,t > 0 , 

d<j>. 


dy 


( t , x, 0) = -w t (t, x ) 0<x<.6,t>0, 


PbWtt + 


d 2 

dx 2 

d 2 

dx 2 


' d 2 w d 3 w \ 

El d ^ + c ° I lh 2 dt ) 


EI~^u{t) [H(x - a n ) - H(x - a l2 )} j 


- p f (f>t(t,x,0) + f(t,x) 0 < x < .6 , t > 0 , 

O r\ 

w (t,0)=—(t,0) = w(t,.6) = -£(t,.6) = 0 t> 0 , 

<^(0 ,x,y) = <j) t (0,x,y) = u>(0,x) = u» t (0,x) = 0 


(5.1) 


where 

f(t,x ) = 2.04sin(1507r/) . 

The parameter choices a = .6 m, i = 1 m, p f = 1.21 kg/m 3 , c 2 = 117649 m 2 /sec 2 , 
p b = 1.35 kg/m, El = 73.96 Nm 2 , c D I = .001 kg m 3 / sec, I< B = 82.9629, T = .0005 m, 
k = 1.9 x 10 -10 m/V, an = .25 and ar ,2 = .35 are physically reasonable for a .6 m by 1 m 
cavity in which the bounding end beam has a centered piezoceramic patch covering 1/6 of 
its length (see Figure 3). The beam is assumed to have width and thickness ,1m and .005 m, 
respectively. The quadratic cost functional parameters were taken to be d\ = d 2 — = 1, 

d 3 = 10 4 and R = 10 -6 with d 3 of much larger magnitude than di,d 2 or d 4 to emphasize 
the penalization of large pressure variations. Note that because there is only one patch, the 
control w T eight R is simply a positive scalar. 

For a beam with the above dimensions and density, the natural frequency of the first 
mode is 73.21 hertz and the frequency of the forcing function was chosen so as to be close 
to this value. To obtain the magnitude 2.04, it was assumed that the forcing function was 
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the result of an exterior plane wave with a sound pressure level of 120 dB (which forces an 
interior sound pressure level of 98 dB). 



0 .6 


Figure 3. Example acoustic chamber with one piezoceramic patch. 


Example 5.1 

In this example, the uniform preservation of exponential stability for the open and closed 
loop approximating systems is examined. For n = m x = m y = 5,6,7 and 8, the margins of 
stability for the open and closed loop systems obtained with tensored Legendre polynomials 
and tensored linear splines are listed in Tables 1 and 2, respectively. The gains needed 
for the closed loop system were calculated via Potter’s method (see [17]). For each n, the 
locations of the open and closed loop eigenvalues obtained with the Legendre polynomials 
are displayed in figures 4 and 5, respectively. When plotting the eigenvalues of A N and 
a n - bvr - 1 (b n ) n N , those eigenvalues having real parts with magnitude greater than 
1 have been excluded in order to better see the distribution near the imaginary axis. Note 
that a uniform margin of stability is maintained between both the open and closed loop 
eigenvalues and the imaginary axis for both sets of bases. Results similar to those obtained 
with the Legendre basis were obtained when finite elements were used as a basis for H™. 
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Table 1 . Margin between the open and closed loop eigenvalues and the imaginary axis with 
tensored Legendre polynomials. 


m x = m y 

n 

max ji?eA,A £ 

max jfleA, A € <r (A" - B N R~ l (b n ) T 11^) J 

5 

5 

-.0145 

-.0196 

6 

6 

-.0213 

-.0220 

7 

7 

-.0200 

-.0200 

8 

8 

-.0158 

-.0290 


Table 2. Margin between the open and closed loop eigenvalues and the imaginary axis with 
tensored linear splines. 














Figure 5. Eigenvalues of 


a n -b n r~ 1 [b n ) t n N 


tensored Legendre polynomials. 


for n = m x = m y = 5,6,7 and 8 with 


As seen in Example 5.1, a larger margin of stability is maintained in both the open and 
closed loop systems with the linear spline basis than with the Legendre basis; hence one might 
conclude that the linear splines are the basis of choice when solving the control problem. As 
noted earlier however, one must also weigh factors such as system size, accuracy and efficiency 
when choosing a numerical method. Numerical tests have indicated that in spite of the larger 
eigenvalue margins of the linear splines, their performance when used in the control problem 
is nearly identical to that obtained with the Legendre polynomials. Moreover, because of the 
lower order accuracy of the splines, a larger number of basis functions is needed to obtain 
suitable accuracy thus leading to matrix dimensions that are almost twice those resulting 
from the Legendre discretization. Finally, as noted earlier in this section, the matrices 
obtained with the Legendre discretization are much more structured than those obtained 
with finite elements or linear splines hence making Legendre implementation more efficient 
than the other cases. Results for the LQR control problem for (5.1) with the tensored 
Legendre basis for the cavity and cubic splines for the beam with m x = m y = 4 and n = 8 
are reported in Example 5.2. 
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Example 5.2 

In this example, the effect of the feedback control on the problem for system (5.1) is de- 
scribed. In order to solve for the optimal control and trajectory, it is necessary to solve both 
the trajectory equation (4.5) and the tracking equation (4.4). Because numerical evidence 
indicated that both unconstrained solutions were roughly periodic with period r = 1/75, the 
problems were solved as initial value problems with starting values y(0) = 0 and r(10/75) = 0 
rather than as free boundary value problems. The choice for initial state is physically rea- 
sonable while the choice to integrate backwards in time in (4.4) is made to reduce numerical 
instability when solving the ODE system for r N (t). 

The uncontrolled and controlled approximate acoustic pressures (p N — Pj<ff) at the 
point (X,Y) = (.3, .1) are plotted in Figure 6 for the time interval [0,10/75]. Similar plots 
for the approximate beam displacement at X = .3 are given in Figure 7. The uncontrolled 
solutions exhibit a beat phenomenon which results from the fact that the frequency of the 
forcing function is slightly greater than the natural frequency of the first mode of the beam. 
After a transient interval, the controlled solutions are periodic and are maintained at a 
level which is approximately 10% of that found in the uncontrolled case (note the scales in 
Figures 6 and 7). This produces an interior sound pressure level of 77 dB which is a 21 dB 
reduction. To further illustrate the state reduction with feedback control, the uncontrolled 
and controlled acoustic pressures at the times T = 1/75, 2/75, 6/75 and 10/75 are plotted 
in Figures 8-11, respectively. The two dimensional plots in each figure show spatial slices 
of the uncontrolled and controlled pressures at X = .3, 0 < y < 1. Figures 12 and 13 
contain plots of the uncontrolled and controlled beam displacements at the times T = 6/75 
and T = 10/75, respectively. The results in Figures 8-13 are representative of those found 
throughout the time interval (0, 10/75] and in conjunction with Figures 6 and 7, demonstrate 
that the pressure and beam displacement are uniformly reduced and maintained at a very 
low level of magnitude in spite of the periodic forcing function. 

The controlling voltage u(t) is plotted in Figure 14. As expected, it is periodic with 
period 1/75. It should be noted that the magnitude of u(i) remains less than 60F which is 
a physically reasonable voltage to put into the piezoceramic patches. 

As mentioned in the last section, the choice of the quadratic cost functional parameters 
di — d 4 and R influences the control stability and performance of the feedback. In this 
problem, the emphasis is on the the reduction of variations in the acoustic pressure; hence 
c ?3 was taken to be larger than d i, d 2 or d 4 . It should be noted that this choice of parameters 
does not exclude the control of the other state variables; in fact, the beam displacement is 
significantly reduced as seen in Figures 7, 12 and 13. Since the parameter R is a penalty 
term for u(t) } more control of the state variables can be effected by choosing R smaller. The 
tradeoff, however, is an increase in the voltage. Hence one must weigh the amount of state 
reduction desired against the amount of voltage which can be put into the patches. 

The amount of control is also directly influenced by patch size, placement and the number 
of patches being used. In the examples that we have observed, the best results were obtained 
with one centered patch, and we have noticed that the amount of control obtained increases 
with increasing patch length. Thus one must weigh the amount of control desired against 
physical limitations on the size of the piezoceramic patches being used. 
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Figure 6. Uncontrolled and controlled pressures at the point (X,F) - (.3,.l) throughout 
the time interval [0, 10/75]. 



Figure 7. Uncontrolled and controlled beam displacements at the point X — .3 throughout 
the time interval [0, 10/75]. 
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xlO 5 T = 6/75 



Figure 12. Uncontrolled and controlled beam displacements at T — 6/75 
uncontrolled displacement, controlled displacement. 

xlO-* T= 10/75 



Figure 13. Uncontrolled and controlled beam displacements at T = 10/75 
uncontrolled displacement, controlled displacement. 
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6 Conclusion 


For the 2-D acoustic problem involving the transmission of exterior noise into an interior 
cavity via fluid/structure interactions, a model set of differential equations has been devel- 
oped. Control is implemented in the model via piezoceramic patches on the beam which are 
excited in a manner so as to produce pure bending moments. By writing the resulting system 
as an abstract Cauchy equation, the problem of reducing interior pressure fluctuations can 
be posed in the context of an LQR time domain state space formulation and approximation 
schemes which are suitable for this theory are presented. 

For one typical patch configuration, examples are given which demonstrate the stabiliz- 
ability of the open and closed loop systems under approximation as well as the reduction of 
cavity pressure and beam displacement when the feedback control is invoked. The examples 
show that input of the optimally controlling voltage u(t ) uniformly reduces both the pressure 
and the beam displacement and maintains them at a very low level of magnitude throughout 
the time interval of interest. 

As mentioned in the example section, the amount of control obtained is directly influenced 
by patch size, placement and the number of patches being used. Initial results have indicated 
that for a uniform periodic forcing function, the best results can be obtained with one 
centered patch with the amount of control increasing with increasing patch length. Further 
computational studies are currently being conducted to determine the effect of various patch 
configurations and forcing functions on the decibel reduction in the cavity. 
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